Colour Fields Computed in SU(3) Lattice QCD for the Static Tetraquark System 
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The colour fields created by the static tetraquark system are computed in quenched SU(3) lattice 
QCD, in a 24^ x 48 lattice at /3 = 6.2 corresponding to a lattice spacing a — 0.07261(85) fm. We find 
that the tetraquark colour fields are well described by a double-Y, or butterfly, shaped fiux tube. 
The two flux tube junction points are compatible with Fermat points minimizing the total fiux tube 
length. We also compare the diquark-diantiquark central fiux tube profile in the tetraquark with 
the quark-antiquark fundamental flux tube proflle in the meson, and they match, thus showing that 
the tetraquark flux tubes are composed of fundamental flux tubes. 
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Figure 1: In the tetraquark flux tube model, the 
elementary flux tubes meet in two Fermat points, at an 
angle of a = 120° to form a double-Y flux tube, except 
when this is impossible and the flux tube is X-shaped. 



Multiquark exotic hadrons like the tetraquark and the 
pentaquark, different from the the ordinary mesons and 
baryons, have been studied and searched for many years. 
The tetraquark was initially proposed by Jaffe [1] as a 
bound state formed by two quarks and two antiquarks. 
Presently several observed resonances are tetraquark can- 
didates. Very recently the Belle Collaboration made the 
tantalizing observation [2], in flve different T(5S) decay 
channels of two new charged bottomonium resonances Zf, 
with masses of lOGlOMeV/c^ and lOGSOMeV/c^ and nar- 
row widths of the order or 15 MeV. Since all standard bot- 
tomonia are neutrally charged, these two new resonances 
have a flavour only compatible with bbud tetraquarks. 
This is the clearest tetraquark candidate so far observed. 
Other potential tetraquark candidates have also been ob- 
served, however they may still be interpreted differently. 
For instance, in 2003, the X(3872) observed by the Bell 
Collaboration [3, 4] was suggested as a tetraquark candi- 
date by Maiani et al [5]. In 2004, the Dsj(2632) state seen 
in Fermilab's SELEX [6, 7] was suggested as a possible 
tetraquark candidate. In 2009, Fermilab announced the 
discovery of Y(4140), which may also be a tetraquark [8]. 
There are as well indications that the Y(4660) could be 
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Figure 2: Tetraquark Wilson loop as defined by 
Alexandrou et al [17], and by Okiharu et al [18]. 



a tetraquark state ]9]. The T(5S) bottomonium has also 
been recently suggested to be a tetraquark resonance ]10]. 
However a better understanding of tetraquarks is neces- 
sary to confirm or disprove the X, Y, Z and possibly also 
other light resonances candidates as tetraquark states. 

On the theoretical side, the first efforts have been to 
search for bound states below the strong decay threshold 
]11-14], as it is apparent that the absence of a poten- 
tial barrier may produce a large decay width to any open 
channel. Recent investigations found that, even above 
the strong decay threshold, the presence of a centrifu- 
gal barrier in high angular momentum multiquarks may 
increase the stability of the system [15, 16[. 

In the last years, the static tetraquark potential 
has been studied in Lattice QCD computations [17- 
19[. The authors concluded that when the quark-quark 
are well separated from the antiquark-antiquark, the 
tetraquark potential is consistent with One Gluon Ex- 
change Coulomb potentials plus a four-body confining 
potential, suggesting the formation of a double-Y fiux 
tube, as in Fig. 1, composed of five linear fundamental 
fiux tubes meeting in two Fermat points [14, 20, 21[. A 
Fermat, or Steiner, point is defined as a junction mini- 
mizing the total length of strings, where linear individual 
strings join at 120° angles. When a quark approaches 
an antiquark, the minimum potential changes to a sum 
of two quark-antiquark potentials, which indicates a two 
meson state. In principle a X-shaped fiux-tube as in Fig. 
lb could also occur, but the potential minimization al- 
ways leads in that case to a two-meson potential. This is 
consistent with the triple fiip-fiop potential, minimizing 
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Figure 3: (a) Lagrangian density 3D plot for ri =8, r2 ~ 14. (b) We also show the 3D plot for ri — 8 and r2 = 8, to 

illustrate that even at distances where the meson meson dominates the flip-flop potential, the meson meson mixing 

with the tetraquark is sufficiently small to produce such a clear a tetraquark double- Y flux tube. The results are 

presented in lattice spacing units (colour online). 



the length, with either tetraquark flux tubes or meson- 
meson flux tubes, of thin flux tubes connecting the dif- 
ferent quarks or antiquarks [14, 16]. 

Here we study the colour fields for the static tetraquark 
system, with the aim of observing the tetraquark flux 
tubes suggested by these static potential computations. 
The study of the colour fields in a tetraquark is important 
to discriminate between different multi-quark Hamilto- 
nian models, quark models with two-body interactions 
only [22], from fiip-fiop models with a multi-body po- 
tential ]16]. Unlike the colour fields of simpler few-body 
systems, say mesons, baryons and hybrids, ]23-26], the 
tetraquark fields have not been previously studied in lat- 
tice QCD. 

To impose a static tetraquark, we utilize the respec- 
tive Wilson loop ]17, 18] of Fig. 2, given by W4Q — 
|Tr (Ml i?i2Af2ii2), where 



TDaa' ^abc^a'b'c' r>bb' tdcc' 

J- aa' abc a'b'c' rbb' j c.c' 



(1) 



The chromoelectric and chromomagnetic fields on the 
lattice are given by the Wilson loop and plaquette expec- 
tation values, 



/p2. ^^ ,p. ^ ^ (W^(ri,r2,r)P(r)o,) 

{E, (r)) = (P(r)„.) - ^w(r,,r,,T)) ^^^ 

{W{n,r„T)Pir),,) _ 

where the jk indices of the plaquette complement the 
index i of the magnetic field, and where the plaquette at 



position r = (x, y, z) is computed at t ~ T/2, 

P^. (r) = 1 - ^ReTr [C/^(r)C/,(r + M)C/t(r + u)Ut{r)] . 

(3) 
The energy (T-L) and lagrangian (£) densities are then 
computed from the fields, 

(H(r)) = l((E2(r)) + (B2(r))) , (4) 

{C{r)) = \ ((E2(r)) - (B2(r))) . (5) 

To compute the static field expectation value, we plot 
the expectation value (^Ef{r)) or (^Bf{r)') as a function of 
the temporal extent T of the Wilson loop. At sufficiently 
large T, the groundstate corresponding to the studied 
quantum numbers dominates, and the expectation value 
tends to a horizontal plateau. In order to improve the sig- 
nal to noise ratio of the Wilson loop, we use 50 iterations 
of APE Smearing with w = 0.2 (as in [25[) in the spa- 
tial directions and one iteration of hypercubic blocking 
(HYP) in the temporal direction. [27[, with ai = 0.75, 
a2 = 0.6 and a^ = 0.3. Note that these two procedures 
are only applied to the Wilson Loop, not to the plaque- 
tte. To compute the fields, we fit the horizontal plateaux 
obtained for each point r determined by the plaquette 
position, but we consider z — for simplicity. For the 
distances ri and r2 considered, we find in the range of 
T G [3, 12] in lattice units, horizontal plateaux with a x^ 
/dof S [0.3, 2.0]. We finally compute the error bars of the 
fields with the jackknife method. 

To produce the results presented in this work , we 
use 1121 quenched configurations in a 24^^ x 48 lattice 
at /3 = 6.2. To test whether these configurations are 
already close to the continuum limit, we first compare 
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Figure 4: Lagrangian density for r2 — 14 and ri from to 6. The black dot points correspond to the Fermat points. 

The results are presented in lattice spacing units (colour online). 
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Figure 5: Colour fields, energy density and Lagrangian density for ri = 8 and r2 = 14. The black dot points 
correspond to the Fermat points. The results are presented in lattice spacing units (colour online). 



the quark-antiquark static potential obtained using these 
configurations with the potential of 381 configurations in 
a larger, 32^ x 64 lattice, at the same f3. The resulting 
quark-antiquark static potentials are identical within the 
statistical error, showing that the volume size effects are 
sufficiently small in our 24^ x 48 lattice. We present our 
results in lattice spacing units of a, with a = 0.07261(85) 
fm or a^^ = 2718 ± 32 MeV. We generate our config- 
urations in NVIDIA CPUs of the FERMI series (480, 
580 and Tesla 2070) with a SU(3) CUBA code upgraded 
from our SU(2) combination of Cabibbo-Marinari pseu- 
doheatbath and over-relaxation algorithm [28, 29]. Our 
SU(3) updates involve three SU(2) subgroups, we work 
with 9 complex numbers, and we reunitarize the matrix. 
We have two options to save the configurations, either in 
a structure of arrays where each array lists a given com- 
plex component for all the lattice sites, or in an array of 
structures where each structure is a SU(3) matrix. 

In our simulations, the quarks are fixed at 
(±ri/2,-r2/2,0) and the antiquarks at (±ri/2,r2/2,0), 
with ri extending up to 8 lattice spacing units and r2 ex- 
tended up to 14 lattice spacing units, in order to include 
the relevant cases where r2 > V^ri. Notice that in the 
string picture, at the line r2 = \/3ri in our (ri, r2) pa- 



rameter space, the transition between the double-Y, or 
butterfly, tetraquark geometry in Fig. la to the meson- 
meson geometry should occur. The results are presented 
only for the xy plane since the quarks are in this plane 
and the results with z j^ are less interesting for this 
study. The flux tube fields can be seen in Fig. 3, 4 and 
5. Theses figures exhibit clearly tetraquark double-Y, or 
butterfly, shaped flux tubes. The flux tubes have a flnite 
width, and are not inflnitely thin as in the string models 
inspiring the Fermat points and the triple flip-flop po- 
tential, but nevertheless the junctions are close to the 
Fermat points, thus justifying the use of string models 
for the quark confinement in constituent quark models. 

In Fig. 6, we plot the chromoelectric field along the 
central fiux tube, (Ey) at a; = 0, for ri = 8, r2 = 14. As 
expected, the chromoelectric field along y is in agreement 
with the position of the Fermat points. The chromoelec- 
tric field along the x — central axis is maximal close to 
the Fermat points situated at a; ~ —4.69 and at a; ~ 4.69, 
flattens in the middle of the flux tube. Outside the flux 
tube, the chromoelectric field is almost residual. 

In Fig. 7, we compare the chromoelectric field for 
the tetraquark and the quark-antiquark system in the 
middle of the fiux tube between the (di)quark and the 
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Figure 6: (-By) in the central axis x = for 

ri = 8, r2 = 14. We show with vertical dashed lines the 

location of the two Ferniat points. The results are 

presented in lattice spacing units (colour online). 



Figure 7: Profile cut at y = of the chromoelectric field 

for the tetraquark and quark-antiquark systems in the 

middle of the flux tube. The results are presented in 

lattice spacing units (colour online). 



(di)antiquark. As can be seen, for our larger distance 
r2 = 14 where the source effects are small, the chromo- 
electric field is identical up to the error bars, and this 
confirms that the tetraquark fiux tube is composed of a 
set of fundamental flux tubes with Fermat junctions. 

To check which of the colour structures, tetraquark 
or meson-meson, produces the groundstate flux tube, we 
study the x^/dof of the T plateaux. Clearly, as expected, 
the X-shaped geometry of Fig. lb never produces ac- 
ceptable plateaux in the range where the meson-meson 
plateaux are good. But, surprisingly, event at distances 
as small as r2 ~ ^tia/S, illustrated in Fig. 3b, where 
the flip-flop potential favours the two-meson flux tube, 
we flnd T plateaux with a good y^ /dof. This shows 
that the mixing between the tetraquark flux tube and 
the meson-meson flux tube is small, and it is possible to 
study clear tetraquark flux tubes even at relatively small 
quark-antiquark distances. 
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